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ABSTRACT 

We present the results of time-resolved optical photometry and spectroscopy of the X- 
ray transient XTE J1859+226 (V406 Vul) . Photometric observations taken during 2000 
and 2008 reveals the presence of the secondary star's ellipsoidal modulation. Further 
photometry obtained in 2010 shows the system ~lmag brighter than its quiescence 
level and the ellipsoidal modulation diluted by strong flaring activity. Spectroscopic 
data obtained with the 10.4-ni GTC in 2010 reveals radial velocity variations of ~ 
500kms~^ over 3h. A simultaneous fit to the photometry and spectroscopy using 
sinusoids to represent the secondary star's ellipsoidal and radial velocity variations, 
yields an orbital period of 6.58 ± 0.05 h and a secondary star's radial velocity semi- 
amplitude oi K2 — 541±70kms~"'^. The implied mass function is f{M) = 4.5±0.6 Mq, 
significantly lower than previously reported but consistent with the presence of a black 
hole in XTE J1859-I-226. The lack of eclipses sets an upper limit to the inclination of 
70 degrees which yields a lower limit to the black hole mass of 5.42 Mq. 
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1 INTRODUCTION 

The mass distribution of compact objects has crucial impact 
on fundamental physics, and can only be determined from 
the s tudy of X-ray binaries (jCharles fc Cod l2006l : ICasarej 
[20o3). In particular, the theoretical black-hole (BH) and 
neutron star (NS) mass distributions depend critically on 
the equation of state (which establishes the critical mass 
dividing NS and BH formation) and on our current unders- 
tanding of the late stages in the evolution of massive stars 
and supernovae. The latter contains several theoretical un- 
certainties which completely dominate the final mass distri- 
bution (e.g. treatment of convective mixing, amount of fall- 
back, mass-cut, mass- loss during Wolf-R ayet phase, initial 
mass function of progenitors, etc; see e.g. iFrver fc Kalogera 
I2OOII ). It should be noted that the upper mas s cut-off of the 
BH d istribution is particularly constraining l|Orosz fc et al.l 
[2003). 
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Only 16 systems out of an estimated Galactic popula- 
tion of ~ 10'^ transient X-ray bi naries have re liable dyna- 
mical mass determinations (e.g. ICasared I2OIC1I ). The typi- 
cal error bars in BH mass measurements are ~40per cent 
and it is clear that more BH discoveries and a factor 2-3 
improvement is neede d before any constraints can be se t 
on supernovae models (|Bailvn et al .111 9981 : [oTel et aLlbOloD . 
The mass determination in X-ray binaries requires knowl- 
edge of the radial velocity curve, the binary mass ratio and 
the orbital inclination (the latter from modelling the tidally 
distorted companion's light curve). 

XTE_J1859+226 was discovered during its 1999 out- 
burst IWood et al. I [1999) and its X-ray properties promptly 
classified it as a BH candidat e. Orbital periods of 6.72 h 
llUemura et al.lll999t). 9.12 h llGarnavich et all 1 19991 ) and 
18.72 h (|McClintock et all I2OO0I ) were reported although 
none of the m could be confirmed . Finally, data taken in 
quiescence (|Zurita fc et al..l |2002| ) during September and 
November 2000 (in 4 night snapshots of 1.5h-4h duration) 
exhibited a ~0.2mag semi-amplitude sinusoidal modula- 
tion, consistent with an secondary star's ellipsoidal variation 
(caused by the changing visibility of the tidally distorted 
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companion star as it orbits the compact object). The peri- 
odogram of this modulation suggests that periods from 6.6 
to 11.2 h are equally pos sible at the 68 per cent confidence 
level. On the other hand, iFilippenko fc Chornockl (fioOli 'l re- 
port a preliminary analysis of radial velocities based on 10 
spectra obtained over two nights. They claim evidence for 
a sinusoidal variation with a period of 9.1 h and companion 
radial velocity semi- amplitude K2 = 570 ± 27kms"\ The 
implied mass function is f{M) = 7.4 ± 1.1 M0, one of the 
largest ever measured making J1859-I-226 a very promising 
massive BH. Although this result has been used by other 
authors, it should be noted that it has never been published 
in a refereed journal and hence we believe these parame- 
ters require confirmation. Therefore, we have embarked on 
a photometric and spectroscopic campaign with the aim of 
determining the orbital period of the binary and its dyna- 
mical mass function. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Photometry 

XTE J1859-f226 was observed on 2008 July 31- August 1 and 
2010 July 13-14 using ALFOSC at the 2.5-m Nordic Op- 
tical Telescope (NOT) in the ii-band. The conditions for 
the 2008 campaign were affected by thick dust (calima) and 
therefore we decided to use long integration times. Eight 
1800-s images were obtained on the night of July 31 and thir- 
teen 900-s images on August 1, covering ~5-7h per night. 
The 2010 campaign was performed under excellent seeing 
(~0.5 arcsec) and transparency conditions, so shorter expo- 
sure times were adopted. Twenty two and thirty two 800-s 
integrations were obtained on the nights of July 13-14 res- 
pectively, covering ~8h per night. XTE J1859+226 was also 
observed on the night of 2010 August 8 using ACAM at 
the 4.2-m Wilham Herschel Telescope (WHT). A total of 
147 frames were obtained in the _R-band using 120-s inte- 
gration time for ~6.5 h. All the images were corrected for 
bias and flat-fielded with IRAF in the standard way. Aper- 
ture photometry was then perfor med on the obje c t and the 
five comparison stars reported in lZurita fc et aU l|2002l ). 



2.2 Spectroscopy 

Optical spectroscopy of XTE J1859-f 226 was obtained on 
the nights of 2010 July 17 and August 13 using OSIRIS at 
the 10.4-m Gran Telescopio Canarias (GTC). Five 19QQ-S 
spectra were obtained every night using the RIOOOB grism 
and a 0.6-arcsec slit, which provides a wavelength coverage of 
AA3630-7500 at 255 km s'^ resolution (FWHM). The radial 
velocity standards 61 Cyg A & B were also observed during 
the second run using an identical set up. A Hg-Ar-Ne arc was 
obtained with the telescope in park position to provide the 
wavelength calibration scale. Standard procedures were used 
to de-bias and flat-field the spectra. The 1-dimensional spec- 
tra were extracted using optimal extraction routines which 
maximize the final signal-to-noise ratio. The wavelength ca- 
libration, flux calibration and radial velocity analysis were 
made with MOLL'lfl. The wavelength calibration was obtained 
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Figure 1. Top: the average optical spectrum of XTE J1859-I-226. 
A red continuum with broad double-peaked emission lines of H« 
and H/3 are clearly present. The feature at ~5890A is possibly 
associated to Hel emission. Bottom: Comparison of the Doppler 
corrected average of J1859+226 with the two observed templates. 
The spectra have been rectified to the continuum and ofli'set ver- 
tically for the sake of clarity. 



using a fourth-order polynomial flt to 34 arc lines resulting 
in a dispersion of 2.14 A pix~^ and a rms scatter of 0.07A. 
The instrumental flexure was monitored by measuring the 
position of the sky line OI 5577. 340A which was used to cor- 
rect the individual spectra. The spectra were also calibrated 
in flux using observations of the flux standard GL57-34. 



3 OVERVIEW OF OPTICAL PHOTOMETRY 
AND SPECTROSCOPY 

3.1 Spectroscopy 

The top panel in Fig. [1] shows the flux-calibrated average 
spectrum of XTE J1859-I-226. It shows a red continuum 
with broad double-peaked emission lines at the positions 
of Hof and H/3, characteristic of X-ray transients in qui- 
escence. The Ha line has a typical equivalent width of 
~130A and FWHM~234Q kms'^ The individual spectra 
of XTE J1859-I-226 and the template stars 61 Cyg A and B 
were prepared for the cross-correlation analysis by subtrac- 
ting of a low-order spline fit to the continuum, after masking 
the Ha and H/3 emission lines and atmospheric absorption 
features. The spectra were then rebinned onto a uniform 
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Figure 2. Top & middle panels: Zurita et al. (2002) epoch C 
optical photometry taken in 2000. Bottom: Optical photometry 
taken in 2008 July 31 and August 1 showing the dominance of 
the ellipsoidal modulation. The solid line is the simulated secon- 
dary star's ellipsoidal modulation determined in section|4] Passive 
refers to the "true" quiescence state as defined by Cantrell & Bai- 
lyn (2007). 



Figure 3. The optical photometry of XTE J1859+226 taken on 
2010 July 13-14 using the NOT (top) and August 8 using the 
WHT (bottom). The short 120-s integration time of the WHT 
images enables us to resolve the flickering activity which strongly 
dilutes the secondary star's ellipsoidal modulation. Active refers 
to the flickering state as defined by Cantrell & Bailyn (2007). 



velocity scale of 112kms~^ per pixel and the unit conti- 
nuum was subtracted. Cross-correlation between the tar- 
get spectrum and the K5V template 61 Cyg A was then 
performed in the region 5000-6516 A after masking the at- 
mospheric and interstellar feature at AA6285-6310 and the 
bump at ~5890 A, possibly associated to Hel emission. Ra- 
dial velocities were extracted by fitting parabolic functions 
to the cross-correlation peaks. Clear velocity excursions of a 
few hundred kms~^ are seen on each night. We note that the 
radial velocities are not significantly affected by the choice 
of either a K5V or K7V template. However, in an attempt 
to constrain the spectral type we have performed an opti- 
mal subtraction analysis of the Doppler corrected average of 
J1859-I-226 (using the orbital sol ution reported in section 
with the 2 templates (see e.g. iMarsh. Robinson fc WoodI 
1 1994 ). The oi the residual slightly favours the K5V tem- 
plate, with 928 versus 944 for 742 degrees of freedom. The 
bottom panel of Fig. [T] displays the Doppler corrected ave- 
rage and the 2 templates. Although the spectrum is too noisy 
to make a direct comparison of the metallic absorption lines, 
the relative depth of the broad CaH and TiO bands at 6300, 
6800 and 7200 A also supports the K5V versus the K7V spec- 
trum. 



3.2 Photometry 

The light curves of XTE J1859-f-226 were determined 
through di fferential photo metry using 5 comparison 



stars (see 
/':'•. fc et al. 



Zurita fc et al] |2002| ) 



as local 
quiescent 



standards, 
magnitude 



1 2OO2D report a 
R=22.48±0.07 but we note that the target is brighter in our 
photometric campaigns by ~0.25 mag in 2008 and ~0.8- 
1.0 mag in 2010. The 2008 data clearly shows an ellip- 
soidal modulation (Fig. [2]), caused by the tidal and rota- 
tional distortion of t he companion star (e.g. see Cen X-4 in 
IShahbaz etlll 119931 ). However, the 2010 light curves show 
considerable scatter due to the presence of strong aperiodic 



variability which distorts the underlying ellipsoidal modu- 
lation: the flickering activity is most prominent in the 2010 
August data (Fig.[3l bottom panel) compared with the 2010 
July campaign (Fig. |3l top panel). This is probably due to 
the shorter integration times which enables us to resolve 
the individual flare events. It is also when XTE J1859-f 226 
is brighter, about 1 mag above quiescence, so it is possible 
that the level of aperiodic variability correlates with the tar- 
get brightness. A similar behaviour is seen in the black hole 
binary A06 20-00, where flickering is foun d to dec r ease w ith 
magnitude (ICantrefl et al.ll2008l : icantreU fc et al.]|20ld ). 

The hght curves of quiescent X-ray binaries are 
known to be contaminated by optical flickering with time- 
scale of min utes to hours and a mplitudes in the range 
0.06- 0. 6 mag (jZurita. Casares fc Sh ahbaz 2003, Hvnes et al.l 
I2OO3I : IShahbaz et al l |2005| . b'oTS ). For comparison, the 
flares in XTE J1859-I-226 have a characteristic time-scale of 
~10min and amplitudes up to 0.5 mag. Therefore, for the 
sake of constraining the orbital period (see section |4]) we de- 
cided to restrict our analysis to the 2008 campaign, when 
the target was closer to its quiescent level and showed mi- 
nimum flickering activit y (the passive state as defined by 
ICantreU fc Bailvnl |2007l). These data were also combined 
wi th quiescent ph otometry obtained in 2000 and reported 
in lZurita fc et alT] (|2002 ) . Both datasets were detrended by 
subtracting the nightly mean value. 



4 PERIOD ANALYSIS AND 

DETERMINATION OF ORBITAL 
PARAMETERS 

To determine the orbital period of XTE J1859+226 we com- 
pute a periodogram using a prior knowledge about the 
expected modulations observed in the binary system. For a 
fixed frequency, we simultaneously fit the photometric Hght 
curve (32 and 21 data points from the 2000 and 2008 cam- 
paigns respectively; Nlc) and the radial velocity curve (10 
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data points; Nrv) with a model to represent the photome- 
tric and radial velocity variations. The photometric ellip- 
soidal modulation is equivalent to sinusoid of frequency 2/ 
(where / is the orbital frequency) which lags by 90 degrees 
in phase relative to a sinusoid at frequency /. Although the 
light curves were taken when the system was in the "passive 
quiescent" state, there is the possibility that the accretion 
disc contribution was not the same for the two epochs, and 
so the light curves will have different amplitudes. To ac- 
co unt for the diffe r ent am plitudes we first fit the light curve 
of lZurita fc et alT] (|2002h and then fit a scaled version of this 
to the data taken in our 2008 campaign. The radial velocity 
motions of the secondary star can be modelled with a sinu- 
soid at frequency /. The phasing of the light curve and radial 
velocity curve are related and so the free model parameters 
are the orbital frequency, phase offset, the semi-amplitudes 
of the three sinusoids, the scale factor to account for the 
different disc contamination between the light curves, the 
magnitude offset for the light curve and the systemic velo- 
city for the radial velocity curve. 

Although the Nyquist frequency sets the maximum fre- 
quency we can theoretically search, inspection of the radial 
velocity curve shows a smooth modulation lasting more than 
2h and so any period we expect to detect should be longer 
than this. Given that observations are taken throughout en- 
tire and consecutive nights, the minimum frequency is li- 
mited to 24 h. Hence we limit our frequency search in the 
range 2 to 6 cycles d^^ with a total of 126 frequency steps 
{2xNt', wh ere NT =NLc+Niiv=63 is the total number of 
data points: [Pressi 2002). Given that there are two different 
types of data with different number of data points, to op- 
timise the fitting procedure we assigned relative weights to 
the different data sets. After our initial search of the para- 
meter space, which resulted in a good solution, we scaled 
the uncertainties on each data set (i.e., the light curve and 
the radial velocity curve) so that the total reduced of the 
fit was «1 for each data set separately. After the scaling, the 
fitting procedure were run again to produce the final set of 
parameters. The 99 per cent white noise significance levels 
was estimated using Monte Carlo simulations. We generated 
light curves and radial velocities with exactly the same sam- 
pling and integration times as the real data, added Gaussian 
noise using the errors on the data points. We computed 5000 
simulated light curves and then calculated the 99 per cent 
confidence level at each frequency t aking into acco unt a re- 
alistic number of independent trial (IVaughanlbOoBI ). 

Three minima above the noise level are clearly seen at 
3.65, 5.21 and 2.69 cycles d~^ (see Fig. 21 bottom panel). 
However, the best fit has a of 54.7 (with 56 degrees 
of freedom) and occurs at a frequency of 3.65 cyclesd"^. 
The best fitted radial velocity parameters for this period 
are listed in Table [T] Note that the K2 velocity, combined 
with the orbital period, implies a mass function f{M) = 
K2 P/2-kG = 4.5 ± 0.6 Mq and hence the dynamical confir- 
mation of a black hole. The next best fit at 5.21 cycles d~^ is 
only significant at the 3x 10~* per cent level and the xevusm- 
in g peak at 2.69 cycles d~^ (i.e. t he ~ 9.1 h period reported 
by iFilippenko fc Chornoc5l200ll ) has an even lower signifi- 
cance. The top panel in Fig. [4] shows the photometric light 
curve and radial velocity curve, folded on the 0.274 d period 
and with the best fitted solution sup erimposed. The peak- to- 
peak amplitude of the light curve of jZurita fc et al..ll2002h is 



P K2 To 70 

(d) (kms-i) (HJD-2455000) (kms-l) 



0.274 ± 0.002 541 ± 70 395.689 ± 0.002 



-28 ± 53 



Table 1. List of the best fitted radial velocity parameters ob- 
tained for XTE J1859-I-226 



0.239 mag and the difference between minima is 0.058 mag. 
The scale factor between the amplitudes of the two light 
curves is 0.84±0.20 which implies that within the uncer- 
tainties, the two light curves have similar amplitudes and 
disc contribution. 

The relatively large amplitude of the light curve and 
the possible difference between the minima (where the mi- 
nimum at phase 0.5 is deeper than the minimum at phase 
0.0) places strong constraints on the binary inclination an- 
gle. Using our X-ray binary model which predicts the light 
curve arising from a Roche-lobe filling in an X-ray binary 
(|Shahbaz et al.l |2003| ) we find that systems with inclina- 
tion angles of ~40 degrees have equal minima. Further- 
more, the lack of eclipses gives an upper limit of ~70 de- 
grees, which implies a lower limit in the mass of the black 
hole of 5.42 Mq. Our data are not sufficiently accurate to 
allow a detailed analysis of the inclination angle, but we 
can use the X-ray binary model to estimate the distance 
and inclination angle. We assume a secondary star with an 
effective temperature similar to that in A0620-00 (an X- 
ray transient with a K4V secondary star with an orbital 
period of 7.8h; see iMarsh Robinson fc Woodlll994l ). E(B- 
V)=0.58(|h^ nes et al.ll2002l ) and the orbital parameters de- 
rived in section 131 We find that to reproduce the large am- 
plitude of the fiux level of the observed light curve requires 
a system at 14kpc with an inclination angle of 60 degrees 
and no accretion disc contribution, or a system with an in- 
clination angle of 70 degrees but with a disc contribution of 
28%. We note in passing that the Galactic rotatio n velocity 
at this distance is ~ — 53kms^ (|Clemens|[l985l ) which is 
consistent with the observed systemic velocity listed in Ta- 
ble [H Clearly a more accurate light curve is required before 
we can place tight constraints on the model parameters. 



5 CONCLUSION 

We have presented 7?-band photometry of XTE J 1859-1-226 
obtained in 2008 and 2010. The 2008 light curves are do- 
minated by the classical ellipsoidal modulation. However, 
the target is ~0.6-l mag brighter than its quiescent level 
during 2010 and the secondary star's ellipsoidal modula- 
tion is found to be strongly diluted by the fiickering ac- 
tivity. A periodogram analysis combining the radial 
velocities and the quiescent photometry yields an orbital 
period of 6.58 ± 0.05 h (low er than previous claims by 
iFihppenko fc Chornockl200ll ) and a velocity semi-amplitude 
541 ± 70kms-^ The implied mass fun ction is 4.5 ± 0.6 Mq 
substa ntially lower than reported in iFilippenko fc ChornocM 
(|200ll) . but still consistent with the presence of a black hole. 
More observations are required to further refine the orbital 
parameters of XTE J1859-f226, in particular, the mass func- 
tion. 
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Figure 4. Bottom: The periodogram. The solid line represents the 99 per cent white noise significance levels. The dashed line shows 
the 99 per cent confidence level above the minimum at a frequency of 3.65 cycles d~^ (=0.274 d). Top left and right: Respectively, 
the optical photometric light curve and radial velocity curve phase folded on the best fit period. 1.5 orbital cycles are shown for clarity. 
In the radial velocity curve plot, the filled circles and stars show the data taken on different nights. In the top light curve plot the filled 
circles are the data from 2008, whereas in the bottom light curve plot the stars are the data from Zurita et al. (2002). The solid line 
in the plots show the best fit model using a double and single sinusoid to simulate the secondary star's ellipsoidal and radial velocity 
variations (see section |4}. 
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